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We present a model for describing flexible polyelectrolytes in a good solvent and in the presence 
of monovalent salt . The molecule composed by monomers is characterized by the end to end 
distance Re = h{Z — l)"*" and the number of associated counterions n. At high temperatures the 
polyelectrolyte behaves as a neutral polymer (7 = 0.588). Decreasing the temperature, the macro- 
molecule changes from this extended conflguration(7 = 0.588) to a stretched form (7 « 1). At 
even lower temperatures, above the Manning condensation threshold, the polyelectrolyte collapses 
(7 ~ 0.3). Our results show good agreement with simulations. 

PACS numbers: Valid PACS appear here 
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I. INTRODUCTION 

Polyelectrolytes represent a very interesting class of 
materials. Biological systems abound with polyelec- 
trolytes. Two of the most well known of all polymers, 
DNA and RNA, are charged polymers QQ- Besides, 
a large class of synthetic polyelectrolytes are present in 
the chemical industry. For instance, polyacrylic acid is 
the main ingredient for diapers Q and dispersions of 
copolymers of acrylamide or methacrylamide and acrylic 
or methacrylic acid are fundamental for cleaning water 
Even though the tremendous interest in polyelec- 
trolyte, unlike neutral polymers 01 @j the understanding 
of the behavior of electrically charged macromolecules is 
still rather poor. 

The contrast between our understanding of charged 
and neutral polymers results from the long range na- 
ture of the electrostatic interactions that introduces new 
length and time scales that render the analytical descrip- 
tion very complicated 0] ■ The presence of charges intro- 
duces more than one new length scale making the scaling 
theories used for no longer as simple and applicable. Be- 
sides the length scale associated with the strength of the 
Coulomb interaction, there is also lengths associated with 
the mean separation of the counterions required by the 
charge neutrality. The association of these free ions to 
the polyion, renormalizing its charge, is another factor 
that has to be taken into account. 

Therefore, simulations seems to be a good technique to 
overcome these major difficulties. Recent molecular sim- 
ulations of salt-free systems 0-0] were able to obtain 
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the end to end distance of strong electrolytes below Man- 
ning parameter [llj . The picture provided by these sim- 
ulations have shown to be more complicated than early 
analytic theories have predicted ^3|"|23- Past theoreti- 
cal works tended to neglect entropy. For stiff chains, such 
as DNA, entropy is an small contribution to the free en- 
ergy and, in principle, can be disregarded. In contrast, 
for flexible polyelectrolytes, treating entropy along with 
with the Coulomb interactions is essential. Acknowledg- 
ing that, one of us developed a theory for describing the 
thermodynamics of salt-free polyelectrolytes in a good 
solvent. Within this approach both entropy and electro- 
static interactions are taken into account 23]^]. The 
end to end distance and pressure calculated usiiig this 
method show good agreement with simulations Q Q [l^ ■ 

The addition of salt is a key factor both in biologi- 
cal and industrial applications. However, including new 
charged species in the system, adds another scale to 
this already very complicated problem. In this work we 
develop a simple theory for studying a dilute polyelec- 
trolyte solution in the presence of monovalent salt in the 
framework of three theories: Debye and Bjerrum for elec- 
trolytes H^-ll^, Manning condensation Tl'l and Flory 
elasticity for polymers 

Unlike neutral polymers, the polyelectrolyte in salt so- 
lution exhibits three distinct behaviors. At high tem- 
peratures the chain is extended. As the temperature 
is decreased, the electrostatic repulsion between the 
monomers along the polyion becomes relevant and the 
chain stretches. At even lower temperatures, above the 
Manning condensation threshold, counterions condense 
into the polyion. The polymer chain net charge decreases 
and the polyion collapses. The addition of salt screen the 
electrostatic interactions and the polymer contracts. Our 
results show a good quantitative agreement with Simula- 
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FIG. 1: The System 

tions Hi]. 

In section II, our model is presented in detail. In sec- 
tion III, the free energy of the system is constructed and 
minimized with respect to two parameters: number of 
ions associated to the polyion and end to end distance. 
The results are presented in section IV. Conclusions end 
this session. 



II. THE MODEL 

We consider a dilute polyelectrolyte solution of con- 
centration pp (see Figure 1). The chains are immersed 
in a continuum solvent with dielectric constant D and 
monovalent salt. There are Z charged groups of diam- 
eter a along the polymer chain, b = 2^l^a is the dis- 
tance between adjacent monomers. The total charge of 
the completely ionized molecule is —Zq < 0, where q is 
the proton charge. The counterions that neutralizes the 
solution have charge +q, diameter a and concentration 
Pm = Zpp. The salt ions have charge +q or —q, diame- 
ter a and density ps- For simplicity the positive ions of 
salt will be refereed as counterions and the negative as 
colons. Thus 

p+ ^ Z{l-m)pp + Ps (1) 
P- = Ps (2) 

are the density of counterions and colons. The inverse 
Debye screening length n is defined in the usual way 
by K = pi Xb, where As = Pq^/D is the Bjerrum 
length, the distance between two ions where the electro- 
static energy equals the thermal energy, pi = p+ -I- p_, 
and (3 = l/kBT. The reduced density for species j is 
p* = Pj<7^, the reduced temperature is T* — Da/f3q^. 
The Manning charge parameter is ^ = j3q^/Db. 
Hence ^ = a/bT* and ^ = Xs/b. For simplicity, 



macroion does not assume any preferential geometry. 
Hence, the distance between monomers i and j is rep- 
resented hy Tij — b\i — jp as usual @. The end to end 
distance Re is the distance between monomers 1 and Z 
and it is represented by: 

Re^biZ~iy. (3) 

The electrostatic interaction between the chain and 
the counterions leads the formation of a complexes made 
of one macromolecule and n associated microions. At a 
given temperature T, monomer density pm = Z/V and 
salt density ps, the system reaches the equilibrium char- 
acterized by a number of associated counterions, m and 
the end to end distance characterized by 7. 

III. THE HELMHOLTZ FREE ENERGY 

The system is, therefore, composed by complexes made 
of one polyion and n counterions, free counterions and 
free colons. The end to end distance of the complex is Re- 
The equilibrium configuration is found by minimizing the 
Helmholtz free energy density with respect to n and Re- 
The free energy density, / = F/Np is composed by three 
contributions, namely 

/ ~ Jelec + fnc + Ient (4) 

where Jelec contains all the electrostatic interactions, 
fnc has the hard core contribution between the different 
species and Jent has the entropic contributions for the 
free energy density. 

The electrostatic free energy density is splited into: 

fsLEC — fpc + /// + fdi (5) 

where fpc accounts for the electrostatic interaction be- 
tween the polyelectrolyte chain and the free electrolyte, 
/// takes care of the electrostatic interaction between the 
counterions e colons and fdi includes the dipole-ion in- 
teraction between the dipole pairs along the chain and 
both charged monomers and free counterions and colons. 

The hard core free energy density is divided into two 
contributions given by: 

fnc = fhc + fcs (6) 

where fhc takes care of hard core interaction between 
monomers along the chain and fcs includes the hard core 
interaction between free ions. Here cs represents the free 
energy of Carnahan-Starling, the hard core term between 
free ions. 

Finally, the free energy density terms that contains 
entropic contributions are: 

fENT = fel + fig (7) 

where fei is the elastic free energy of the flexible chain, 
fig is the ideal gas free energy for the mixture of the 
different species present in the solution. 
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A. The polyion-free ions electrostatic free energy 
density fpc 



the usual dipole-ion interaction originally derived for the 
coulomb gas and given by [s^, |^ : 



The electrostatic free energy density between the com- 
plex made of a polyelectrolyte and n counterions associ- 
ated to it and the free ions in the solution can be derived 
in the framework of the Debye-Huckel theory yielding ( 
see appendix A): 



Pipe 



where Pz — ~^ 



m — — 1 + Z/n and 

-Kr{x) 



1=2 1 dx{Z -x) 

10 



r{x) 



(8) 



(9) 



where r{x) = b x'^ . 

These ideas have been successful in describing charged 
systems, including simple electrolyte solutions l30l ISl f . 
charged rods :32, 33j, char g e d sp herical colloids j3J|, and 
flexible charged chains [23, |2J, y^- In particular, the 
binding isotherms of DNA with surfactant molecules have 
been found in good agreement with experiment [3^ . 

The integral X might also be writen using the incom- 
plete gamma function 'ST] , namely 



7mc(a,a;) 



dt e~H' 



t,a~l 



Rea> . 



(10) 



The sufSx inc is used to avoid confusion with the expo- 
nent 7. In this case, we get 
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Zk 



K 



7mc(-l + l/7,«^'^^) (11) 

7„c(-l + 2/7,«fc^^) 

^2-7 



6(1-7) 6(2-7) 



(12) 



B. Free ion- free ion electrostatic free energy 
density/// 

For the interaction between free ions we employ the 
Debye-Hiickel free energy for electrolytes 26, 27], 



1 



[ln(l + Ka) - KO- + (Kcr)V2] • (13) 



C. The dipole-ion electrostatic free energy density 

fdi 



1 



Zm{cj/a2Y x^uj2{x) , (14) 



where x = ^02, kg = \j , p*_ = Z (1 — m) cr'^/V'e, 
Ve = 4 7ri?3/3, 02 = 1.1619(7, and 

uj2{x) = 3[ln(l + X + xV3) -x + x^lQ\/x'^ . (15) 

Note that we use here k in order to take the interac- 
tions between the dipoles formed on the chain and the 
charged monomers and free ions. It is essential for the 
attractive forces that lead to collapse of the chain at low 
temperatures. The density pl_ is the density of charged 
monomers not associated in the volume occupied by 
the chain. 



D. The excluded volume free energy densities f^c 

and fcs 

The polyelectrolyte will be consider to be in a good sol- 
vent. Therefore, hard-core repulsion between monomers 
is approximated by a virial coefficient from Flory-de 
Gennes theory H lE H N , 



(16) 



where Wi = Awa^ /3 is the second virial coefficient for 
hard spheres of diameter ct, and p = Z/Ve with Ve = 
47ri?3/3. 

The hard-core repulsion between free ions is approxi- 
mated by the Carnahan-Starling free energy |40|. 



f3fcs = V piy 



4-3;/ 

(1-y)' ' 



(17) 



where y = i^pX/Q is the volume fraction occupied by the 
free positive and negative ions. 



E. The elastic free energy density /ej 

We are considering here flexible polyelectrolytes. The 
entropic contribution of the elastic free energy densit y is 
the same of a neutral polymer and it is given by ja, l35l 



When a counterion associate to a negative monomer 
of the polyelectrolyte, it forms dipole. The interac- 
tion between these dipoles formed by association and 
the monopoles consisting of free counterions and colons 
and non-associated monomers along the chain is given by 



(18) 



where a = Re/Ro, where Rq = h\/ Z — \ is the end to 
end distance of a polymer in a O solvent. 
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F. The ideal gas free energy density fig 



IV. RESULTS AND CONCLUSIONS 



Another important cntropic contribution is due to the 
mixing of the different species, complexes, counterions 
and coions. The free energy density associated with the 
mixing of ideal particles is given by [s^, |^ : 

pf,g = — [p+lnp;A3-p++p_lnp*_A3 (19) 
Pp 

- p_]+lnppA3-l + /5/4 



where A = h/\/2TTmkBT is the thermal wavelength, /f^ 
includes the internal degrees of freedom of the complex 
made of a polyelectrolyte and the associated ions and it 
is given by: 



= Z(l-m)[ln 



Z{1 - m) 



A'^-ll 



(20) 



+ Zm[\n§!^A'-l]-Z[\n^A'-l] 

VeC,2 Ve 



Since the monomers along the chain are mobile, the 
first two brackets in Eq. (|20|l account for the entropy of 
mixing of charged monomers and dipoles along the chain. 
The third bracket discounts the overcounting since the 
entropy of a neutral polymer was already accounted in 
Eq. H18|) . The last term is the internal partition function 
of the complex in its internal reference frame. The in- 
ternal partition function of the dipoles, ^2, is given by: 



C2 



47r 



e' 



'/'■^* dr = K{T) , 



(21) 



where Rsj is chosen according to Bjerrum to be a /2T* , 
the location of the minimum of the integrand. Thus the 
association constant is different from zero for T* < 0.5. 
Above this temperature there is no countcrion association 
and no dipole formation on the chain. We have then 

IMIsSEil 



(22) 



with 



Qo 



-L^e-'/^'mi/T*)-m2) + e^] (23) 
1_ /J_ J_ 

The internal partition function of a n-complex con- 
tains: one electrostatic term due to the interaction be- 
tween the monomers and associated ions and one term 
due to the different ways n ions can associate into Z 
monomers in a polyions. The additions of these two 
parcels gives: 

z-i 

Pf^' = ^plbY.{Z-t)— (24) 



i=i 



+ Z m In m + Z (l — m) hi{l ~ m) 
with r(i) — bi'^. 



The behavior of the polyelectrolyte under the variation 
of temperature is illustrated in Figure 2. At high tem- 
peratures ( low Bjerrum length) the end to end distance 
approaches the one of a neutral polymer in a good sol- 
vent and Re = b{Z - 1)0-588, p^j. z = 32 and b = 2^/^a 
the end to end distance is given by Re w 8.45f7. As the 
temperature decreases, the electrostatic interactions be- 
come dominant, the polymer stretches and the end to 
end distance is given by 7 « 1. Above the Manning 
condensation threshold, the counterions associate to the 
monomers along the chain, decreasing the electrostatic 
interactions. The elastic energy becomes dominant and 
the polymer collapses. Our results are compared with 
simulations (2^. The actually salt density used for that 
simulation is not specified in the manuscript but is within 
the interval p* = 0.001 and p* = 0.008. 

Figures 3 and 4 show that for a given Z , the radius of 
end to end distance, Re depends only on the value of the 
bare inverse of screening length, kq = ■\/47r(pm -I- 2ps)Xb, 
and not on the density of polyelectrolytes. In Figure 3 
the renormalized end to end distance given by 



R* = 



Re 



6(Z/32)0-588) 



(25) 



is ploted against In kqct for Z = 32, 64 for varius densities 
and for Xb = 0.83(7 (^ = 0.75). The decrease of the end 
to end distance is, in this case, due to the screending of 
the salt present in the solution. Note that here there is no 
condensation. The renormalized end to end distance, R* , 
was constructed as follows. R was divided by the end to 
end distance of an equivalent neutral polymer in a good 
solvent Ro « 6^°-^^*. In order to keep both Z = 32 and 
Z = 64 in the same scale, N was also divided hy Z = 32. 
Our results agree with simulations performed in the scale 

In Figure 4 we show R* for As = 3.2cr(^ = 2.9). In this 
case, since the system is above the Manning threshold, 
condensation is present and the polyelectrolyte should 
form a more compact form than for Xb = 0.83(7 (^ — 
0.75). Comparison with Figure 3 shows that for Iukqct > 
—0.5 larger Xb has larger R. Part of this effect is due to 
plotting versus Iukoct. For the same p„i and p^, the end 
to end distance is surely smaller for Xb — 0.83(7 than for 
Xb = 3.2(7(^ = 2.9) as indicated by Figure 2. 

In summary, we studied the end to end distance of 
polyelectrolytes with added salt. We found that R for 
a fixed Z and Xb depends only on the amount of salt 
in solution. In this case, the main effect of the elec- 
trolytes is to screen the interactions and decreasing R. 
For fixed values of Z, pp and ps, as the temperature is 
decreased, R exhibits three distinct behaviors: stretched, 
extended and collapsed. This three regimes are related 
to the stretched polymer, fully extended polyelectrolyte 
and collapsed polymer behavior. 
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FIG. 2: Plot of end to end distance as a function 
of As for Z = 32, = 0.001, and from top 
to bottom = — 0.008. The circles are the 
simulational results extracted from ref. 12911. 
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FIG. 4: Same as in Figure 3 but for ^ = 2.9. 
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In OKq 

FIG. 3; Plot of the renormalized end to end 
distance 7?* as a function of In ko for Z = 32 
( solid line ) and Z = 64 { dashed line), for 
p'^ = 0.0001,0.001,0.01,0.02 and ^ = 0.75. The 
points for different densities fall in the same line. 
The symbols are the simulation results for the 
same pj^ and Z = 32 (cir cles) and Z — 64 
(squares) extracted from ref. |29|. 
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APPENDIX A: CALCULATION OF fp^ 

In order to compute the free energy density associated 
with the polyion-free ions interaction, let us start by cal- 
culating the electrostatic potential at a point P due to 
the charged chains and electrolytes. Each chain has a 
linear charge density A' — —{Z — n)q/L. The coordi- 
nates on the molecule are denote by a prime, and the 
coordinates of the point P where we want to evaluate 
the electrostatic potential are P = {x,y,z). The electro- 
static potential at a point P located in an ionic solution, 
due to a charge element X'ds' on the molecule, is given 
by the Debye-Hiickel expression, 

-Kr 

— , (Al) 



Dr 

where r is the distance of the charge element to the point 
P. Integrating the expression above, the electrostatic 
total potential at P becomes: 



D 



ds'- 



(A2) 



where L = b{Z — 1) is the contour length. The electro- 
static potential at P due to the chain when there is no 
electrolyte is given by: 



A' f 1 
V - - / ds'- , (A3) 



the suffix pr meaning proper. 

The potential due to the ionic solution discounting the 
chain is given hy tjj = (j) — (f'pr- The free energy is then 
evaluated by using the charging process of Debye [2^ 
in this potential. 

The electrostatic energy of a charge element dq at 
(x, y, z) in the potential ip is 



dU = dq{x,y,z)ij{x,y,z) 



(A4) 
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The electrostatic energy of the entire molecule is obtained 
integrating the above expression, what gives 



^ rL 
2D 



where r = r(s — s') = r(s' — s) denotes the distance 
between line elements ds and ds', both on the molecule. 
The factor 1/2 is included to avoid double counting. Now 
we substitute q by <^q in U and integrate in from to 
1. We obtain 



Therefore, 



v2 fL 



(A6) 



(A7) 



L= ds= bdt . 
Jo Jo 



Hence, 



pc 



r{t - 1') 



2D 



f dt [ 

Jo Jo 



dt' 



(A8) 



(A9) 



Changing variables again, x = t — t', and integrating by 
parts, the free energy becomes 



fpc — - ^ ^ 1 1=2 



2D 



[ dx{Z- x) 
Jo 



g-ftr(x) _ ^ 



r{x) 



(AlO) 



Defining pz — —l + m and substituting A = (— 1 
m)q/b = Pzq/b, we may write 



pc 



(All) 



We now change from s to a dimensionless variable t, with 
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